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1 INTRODUCTION 



In the last few years there has been enormous progress in string theory, connected with various 
ideas of duality. Similarly in N = 2 supersymmetric (susy) gauge theory the idea of electromagnetic 
(EM) duality emerged in 1994 with the work of Seiberg-Witten ||30|| . The mathematical consequences 
of such work have been no less profound than the physical import but a survey is not contemplated 
here. One theme which has emerged however involves the interaction of ideas of integrability with 
matrix models, topological field theory (TFT), EM duality, Landau-Ginsburg (LG) models, etc. 
This is perhaps not surprising in that integrability involving invariant tori etc. seems fundamental 
in quantization itself but the interaction here is most explicit and inevitable. Without going into 
a discussion of N = 2 and N = 1 susy let us simply say that the low energy effective action 
for an ./V = 2 susy Yang-Mills (YM) theory with gauge group SU(2) for example is described by a 
holomorphic prepotential T (see ]|, |9|, ^] for a discussion and |19| for susy - we will not give extensive 
references to the physics here). This theory has a scalar potential V(4>) = in the vacuum and there 
are generally nonvanishing <p f° r which this occurs. Setting <p = (l/2)ao"3 (03 = ( J -1 )) an< ^ 
u =< Trcfr 2 >= (l/2)a 2 (< >~ vacuum expectation value or vev, where ~ is used for "corresponds 
to" or sometimes "is asymptotic to"), the complex parameter u labels inequivalent vacua and the 
manifold of gauge inequivalent vacua is called the moduli space M of the theory (here u is a coordinate 
on Ai ~ C with punctures or singularities). A metric on M is given locally by ds 2 = QJ-"(a)dada = 
Qr(a)dada where r(a) ~ (9/2n) + (4-7ri/g 2 ) is the complex coupling constant. One writes a D = 
(dJ- '/da (with < <p D >= (l/2)a D a^) and the superscript D indicates duality here. There is a 
Legendre transformation connecting dual superfields, etc. and ds becomes ds 2 = Qda D dd. On 
M. physical arguments show that a(u), a D {u) are single valued except for monodromies around 

say (uq, —uq, 00) (e.g. (a D ,a) T — > M QO (a D ,a) T for = ( ^ -1 )) an< ^ one P^ks = 1 by 
renormalization arguments. This leads to a realization of the monodromies by a differential equation 
(4&) {—D 2 + V(z))ij}(z) = 0, ip = ((i/2)a, — ia D ) T and this can be rewritten as a hyper geometric 
equation from which one obtains 

,u — l A 1 1 — u V2 f u dx^Jx - u 
a D {u) = iip 2 (u) = i——F{-,-,2,——) = / — • (1.1) 
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a(u) = -2iV»i(u) = v / 2(n + l)2F(--,-,l, ' 
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One can eliminate u to get e.g. a D (a). The integral formulas in (1.1) are easily recognized as arising 
from an elliptic curve (take branch cuts (—1,1) and (it, oo) in C and construct a two sheeted Riemann 
surface of genus one (torus)). Take homology cycles A ~ cycle 1 — > u on one sheet and u — > 1 on the 
other, with B ~ cycle around (—1, 1) and define now 



dX; 



n 



B 



dX; dX 



1 



u 
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zdz 



(1.2) 



This differential dX = dXsw is the Seiberg-Witten (SW) differential and the integrals in (1.2) are 
period integrals which satisfy (<(k) under the label of Picard-Fuchs (PF) equations. The correct 



monodromies are obtained directly (cf. also [32] for an interesting use of monodromy ideas). We will 
indicate how all this is related to integrable systems and the Whitham equations. 



2 WHITHAM EQUATIONS 

We will work in the context of algebro-geometric solutions of integrable partial differential equa- 
tions (PDE) which involves the following ingredients (cf. H ||, 14] ). Take an arbitrary Riemann 



surface £ of genus g, pick a point Q and a local variable 1/k near Q such that k{Q) = oo, and, for 
illustration, take q{k) = kx + k 2 y + k 3 t. Let D = Pi + • • • + P g be a non-special divisor of degree 
g and write ip for the (unique up to a constant multiplier by virtue of the Riemann-Roch theorem) 
Baker- Akhiezer (BA) function characterized by the properties that ip is meromorphic on £ except 
for Q where ip(P)exp(—q(k)) is analytic and (*) ip ~ exp(q(k))[l + Y^TiCj/k-')] near Q; on S/Q, ip 
has poles at the Pi. In fact ip can be taken in the form (P G E, Pq ^ Q) 

t .p) = ^ / . e( - 4(F) ;^:f: + " v,+z ° ) (2.D 

JPo <c>{A{P) + z ) 



where dtt 1 = dk + • • • , dft 2 = d(A: 2 ) + • • • , dQ 3 = d(k 3 ) +---,Uj = f B . dn 1 , Vj = J B dn 2 , Wj = 
J B . dQ 3 (j = 1, • • • ,g), zq = —A(D) — K, and G is the Riemann theta function. Here the dflj are 
meromorphic differentials of second kind normalized via j A ddj = (Aj , Bj are canonical homology 
cycles) and we note that xdQ, 1 +yd£l 2 + tdQ 3 ~ dq(k) normalized; A is the Abel-Jacobi map (A(P) = 
J P diOk where the dujk are normalized holomorphic differentials, k = 1, f A .du>k = 5jk), and 

K = (Kj) ~ Riemann constants {2K = —A{Ky) where K a is the canonical class of £ ~ equivalence 
class of meromorphic differentials) so Q(A(P) + Zq) has exactly g zeros (or vanishes identically). 
The paths of integration are to be the same in computing Jp dQ l or A(P) and it is shown in 
[||, ^, 14] that ip is well defined (i.e. path independent). Then the £j in (*) can be computed 



formally and one determines Lax operators L and A such that d y ip = Lip with dtip = Aip. Indeed, 
given the £j write u = —2d x ^\ with w = 3£,id x £,i — 3<9 2 £i — 3d x £,2- Then formally, near Q, one 
has (-d y + d 2 + u)ip = 0(l/k)exp(q) and (-d t + d 3 + (3/2)ud x + w)tp = 0(l/k)exp(q) (i.e. this 
choice of u, w makes the coefficients of k n exp(q) vanish for n = 0, 1, 2, 3). Now define L = d x + u 
and A = d 3 + (3/2)ud x + w so d y ip = Lip and dtip = Aip. This follows from the uniqueness of BA 
functions with the same essential singularity and pole divisors (Riemann-Roch). Then we have, via 
compatibility L t — A y = [A,L], a KP equation (3/4)u yy = d x [u t — (l/4)(6uu x +u xxx )\ and therefore 
such KP equations arise automatically from a Riemann surface and are parametrized by nonspecial 
divisors or equivalently by points in general position on the Jacobian variety «/(£)■ The flow variables 
x, y, t arise via q(k) and then miraculously reappear in the theta function via xU + yV + tW; thus 
the Riemann surface itself contributes to establish these as linear flow variables on the Jacobian 
and in a certain sense defines the flow variables. The pole positions Pi do not vary with x, y, t and 
(t) u = 2d 2 loge(xU + yV + tW + z ) + c exhibits G as a tau function. 
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We recall also that a divisor D* of degree g is dual to D (relative to Q) if D + D* is the null divisor 
of a meromorphic differential dCl = dk + {j3/k 2 )dk + • • • with a double pole at Q (look at £ = 1/k to 
recognize the double pole). Thus D + D* - 2Q ~ £T S so 4(£>*) - _4(Q) + if = -[.4(D) - .A(Q) + K]. 
One can define then a function tp*(x,y,t, P) = exp(—kx — k 2 y — k 3 t)[l + • • •] based on 

D* (dual BA function) and a differential <if2 with zero divisor D + D*, such that (j) = i/jijj*dCl is 
meromorphic, having for poles only a double pole at Q (the zeros of dCl cancel the poles of ipip*). 
Thus ipip*d(l ~ ipip*(l + ((3/k 2 + ■ ■ -)dk is meromorphic with a second order pole at oo, and no other 
poles. For L* = L and A* = -A + 2w - (3/2)u x one has then (d y + L*)i(j* = and (d t + 4*)V>* = 0. 
Now the prescription above seems to specify for tp* (U = xU + yV + tW, Zq = —A(D*) — K) 

* - J? { X d&+yd&+tdn 3 ) Q(A{P) -U + zg) 

In any event the message here is that for any Riemann surface E one can produce a BA function 
ip with assigned flow variables x,y,t,--- and this ip gives rise to a (nonlinear) KP equation with 
solution u linearized on the Jacobian J(S). 

One knows that hyperelliptic curves play a special role in the theory of algebraic curves and 
Riemann surfaces. For hyperelliptic Riemann surfaces one can pick any 2g + 2 points Xj € P 1 and 
there will be a unique hyperelliptic curve T, g with a 2-fold map / : T, g — > P 1 having branch locus 
B = {Xj}. Since any 3 points Aj, Aj, A& can be sent to 0, 1, oo by an automorphism of P 1 the 
general hyperelliptic surface of genus g can be described by (2g + 2) — 3 = 2g — 1 points on P 1 . 
Since / is unique up to an automorphism of P 1 any hyperelliptic T, g corresponds to only finitely 
many such collections of 2g — 1 points so locally there are 2g — 1 (moduli) parameters. Since the 
moduli space of algebraic curves has dimension 3g — 3 one sees that for g > 3 the generic Riemann 
surface is nonhyperelliptic whereas for g = 2 all Riemann surfaces are hyperelliptic (with 3 moduli). 
For g = 1 we have tori or elliptic curves with one modulus r and g = corresponds to P 1 . In 
many papers on soliton mathematics and integrable systems one takes real distinct branch points 
Aj, 1 < j < 2g + 1, and oo, with Ai < A2 < • • • < X2 g +i < 00 and /u 2 = — ^j) = ^23+1 (A, Aj) 

as the defining equation for E g . Evidently one could choose Ai = 0, A2 = 1 in addition so for g = 1 
we could use < 1 < u < 00 for a familiar parametrization with elliptic integrals, etc. One can take 
dX/ fj,, XdX/ fj,, • • • , X 9 ~ 1 dX/ ' fi as a basis of holomorphic differentials on E g but usually one takes linear 
combinations of these denoted by duj, 1 < j < g, normalized via § A . du)j = 6ij, with period matrix 
defined via § B . duoj = bij. The matrix II = (Hij) is symmetric with 3TI > and it determines the 
curve. One will also want to consider another representation of hyperelliptic curves of genus g via 
fi 2 = IIo 9+1 (^ — -\?) = P2g+2(X, Xj) where 00 is now not a branch point and there are two points n± 
corresponding to A = 00. 



We extract now from |7]] which provides some embellishments to [20, 25]. Variations of the 



Whitham-Bogoliubov averaging procedure have proved to be of considerable interest and a very nice 



introduction to averaging for KdV appears in [18]. Thus one considers modulated waves based on a 



potential u(9i, Xj) where the fast variables x, t occur via 9i = xUi+tWi and the slow variables occur in 
the moduli (here branch points) Aj = Xj(X,T) (for KdV one is dealing with hyperelliptic surfaces). 
The averaging procedure looks for example at < f(u) >~ Zim/,_ >00 (1/2L) J_ L f(u){dx or dt} for 
suitable functions f(u), and, assuming incommensurable frequencies, ergodicity implies < f(u) >= 
(1/2tt) n Jq 77 ■ ■ ■ Jq 77 f(u(0{, Xj) JJi N d6i, where the slow variables are fixed. This is made more explicit 
below. One should perhaps contrast this procedure explicitly with dispersionless limit techniques 
which arise in the genus zero situation. For KdV 114 + (3/2)1^ — (l/4)u xxx = (based e.g. on 
L 2 = d 2 — u, B = d 3 - {3/2)ud - (3/4)u x , d t L 2 = [B, L 2 ]) one writes et = T, ex = X to obtain 
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formally ut + (3/2)uux = as e — > (dispersionless KdV or Euler equation). The background 
dKdV theory for Lax operators involves Lip = \ip with xfi = exp[(l/e)S(X,T)] so for P = dxS, A = 
P 2 — U(X,T) where one usually assumes u(x,t) = U(X,T) + 0(e). This last step is never really 
discussed adequately in work on (algebraic type) dispersionless theory and we make a brief comment 
here. In question is the behavior of u(x, t) = u(X/e, T/e) as x, t — > oo (or e — ► with X, T fixed) and 
there is a priori no modulated wave or quasiperiodic situation to make the analysis easier. A detailed 
analysis of KdV — ► Euler for example is developed in classical work of Lax, Levermore, and Venakides 
for example (see [^] for references and discussion) and a weak limiting procedure is justified in many 
situations. Similar analysis also applies to nonlinear Schrodinger and Toda equations. There is 
however no general theory here and the analysis becomes deep and technical. The moral is that the 
algebraic limiting procedure is justified sometimes in realistic situations and we simply assume it to 
be OK whenever it arises. The averaging procedure for quasiperiodic modulated waves on the other 
hand seems to be quite generally applicable. 

Let us sketch now the averaging procedure following p3| with clarifications as in |7|, |2(| (cf. also 
). Consider KP in the form 3u y y + d x (Aut — Quu x + u xxx ) = via compatibility [d y — L,df — A] = 
where L = d 2 — u and A = d 3 — (3/2)ud + w (u — > — u in (2.1). We have then (d y — L)ip = 



with (dt — A)ip = and for the adjoint or dual wave function 0* one writes 0*7 = —d y ^* with 
ip*A = -d t ip* where 0*(/#O = (-<9) J '(^7)- We can use formulas (cf. Q and Q) = 
e px+Ey+m -(t)(Ux+Vy+Wt,P) and tp* = e - px - Ey - nt -ct>*(-Ux-Vy-Wt, P) to isolate the quantities 
of interest in averaging (here p = p(P), E = E(P), Q = £l(P), etc.) We think here of a general 
Riemann surface T, g with holomorphic differentials duo^ and quasi-momenta and quasi-energies of the 
form dp = dn 1 , dE = dn 2 , dQ = dQ 3 , ■■■ (p = Jp Q dtt 1 etc.) where the dQ j = dttj = d(\ j + 0(\- 1 )) 
are meromorphic differentials of the second kind. Following j25| we normalize now via 3ft J A dVL k = 
5ft f B .dfl k = 0. Then write e.g. U k = (l/2m) § Ak dp and U k+g = — (l/2m) § Bk dp (k = 1, ••-,(/) 
with similar stipulations for V k ~ <fd£l 2 , W k ~ §d£l 3 , etc. This leads to real 2g period vectors 
and evidently one could also normalize via § A dQ k = or 9 6* dQ k = 9 $ R d£l k = 0. It is then 
immediate that (ip*L)ip = ■■ where e.g. ^*(u&>) = (-d) j (ip*u). 

More precisely, for A = Eo a i d \ one has V " A = zZoi-df^*^) and (^*A)ip = £o &> (ip* (A^ ») 
with A = A, A 1 = - YXiaid 1 ' 1 , A 2 = Y%W ~ l)/2]a i d i - 2 ,- ■ ■. For L and A as in KP above 
L 1 = -2d, L 2 = 1 and A 1 = -3d 2 + (3/2)u, A 2 = 3d, A 3 = -1. 

Now for averaging we think of it = u ([l/e]S(X, Y, T)\I(X, Y,T)) + eu x (x, y, t) + e 2 u 2 (x, y, t) H 

(7 ~ moduli and S ~ action in some sense) with dxS = U, dyS = V, and dxS = W. We think of 
expanding about uq with d x — > d x +edx- This step will cover both x and X dependence for subsequent 
averaging. Then look at the compatibility condition (J) : dfL — d y A + [L, A] = 0. We will want the 
term of first order in e upon writing e.g. L = Lq + tL\ + • • • and A = Aq + eA\ + • • • where Lq, Aq 
are to depend on the slow variables X, Y, T, and this gives dtL\ — d y A\ + [Lq, Ai] + [L±, Aq] + F = 
0; F = dxL — dyA — L 1 dxA + A 1 dxL. Thus F is the first order term involving derivatives in the 
slow variables. Next one writes d t (ip* (Li^)) - dy(ip*(Ai-iJj)) = t/j*{L lt - A ly + [Xi, A] + [7,^4i]}V> = 
tp*(dtLi — d y A\ + [LqjAi] + [Li,Aq])i[) + d x (- ■ •) and via ergodicity in x,y, or t flows, averaging of 
derivatives in x,y, or t gives zero so one obtains the Whitham equations in the form < ip*Fijj >= 
(this represents the first order term in e - the slow variables are present in Lq, Aq, i/j, and ■0*). In 
order to spell this out in [^] one imagines X, Y, or T as a parameter £ and considers 7(£), A(£), 
etc. (in their perturbed form) with 0(C) = e P^ x+E ^ y+n ^ t ■ </>(U(£)x + V(£)y + W{Qt\I(£)) and 
ip* = exp(—px — Ey — Qt)cj)*(—Ux — Vy — Wt\I) (no £ variation - i.e. assume p, E, Q, U, V, W, I fixed). 
Also note that x,y,t and X,Y,T can be considered as independent variables. Upon differentiating 
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various expressions with respect to x, y, t and £ and combining one obtains then 



dfl < ip*ip >= -dp < ^A 1 ^ >; dE < ?p*ip >= -dp < ^*L l ij} > (2.3) 

and a version of the Whitham equations in the form 

p T = fl x ; p Y = E x - E T = fly (2.4) 

Note also from dxS = U, dyS = V, and 8tS = W one has compatibility relations dyU = 
dxV, dxU = dxW, and 8tV = dyW. We mention also an equation obtained along the way which 
will be important in the analysis of pq|, namely 



< tp* (L^dxA - A l d x L)ip >= d x fl < V*^V > ~dxE < f > + 

+ d x W- < ip*L l ip e > -d x V- < ip*A l ip e > (2.5) 

We feel that this derivation is important since it exhibits again (as in [|l^] for the KdV situation) the 
role of "square eigenfunctions" (now in the form ip*ip) in dealing with averaging processes (cf. M, |T(| 
for more on ip*ip and [jjj for details on averaging). 

Next, d'apres [Jjj, we give some formulas for differentials followed by a few remarks on hyperelliptic 
situations. The KP flows can be written as d n u = K n (u) where the K n are symmetries satisfying (in 
the notation of fK|) the linearized KP equation d 3 /3 = (l/4)<9 3 /? + 3d(u(3) + (3/4)a" 1 5|/3 = K'[f3]. 
The conserved densities or gradients 7 satisfy the adjoint linearized KP equation ^37 = (l/4)<9 3 7 + 
3ud^/ + (3/4)9 -1 <927. Then, replacing the square eigenfunctions of KdV theory by ipip* one has e.g. 
V'V'* = Z)o° s n^~ n where s n ~ j n . Further d n u = K n+ \ = ds n+ \ = dRes L n = dVI n+ i where 
^/ ~ 3f/3 u (Res L n = nH^-i is generally used in the multipotential theory). We are thinking 
here in a single potential theory where all potentials U{ in L = d + u i+id~ l are expressed 
in terms of 112 = u via operations with d and d . One uses here the Poisson bracket {f,g} = 
I J (Sf /5u)d(5g/5u)dxdy (Gardner bracket). Thus one has s n+ i ~ 7 n +i ~ Vi n +i as conserved 
gradients satisfying the adjoint linear KP equation (**) dtj = (l/4)<9 3 7 + 3u^7 + (3/4)<9 _1 <9^7. 
Following |l8|] one linearizes around a fixed finite zone solution u in the adjoint linear KP equation 
and puts this u into the h etc.; then averaging over the Oi variables is carried out as before. Also from 
il>ip* = X)i°(sn/A™) (cf. 0) we can write < ipip* >= Y,T(< s n > A") Such a series is natural from 
asymptotic expansions but when tp, ip* are written in terms of theta functions it requires expansion 
of the theta functions in 1/A (such expansions are documented in [13], p. 49 for example). It is 
now natural to ask whether one can express dp, dE, dfl from (|273 ) in more detail. Some order 
of magnitude considerations suggest dp ~< ipip* > d\ near 00. To prove this we refer to [1C] 



for background and notation (cf. also JTTj, [TJ] for dispersionless genus zero situations). We write 
8 = L + YT <j)L~ j which implies that (dip/ip) = A + £i° <rjA _J '. Then (logiP) x = X + and 
(logip) x ~< (logip) x >= p = A + J2T < a ] > ^ ^ ■ But s n+ i = —no\ — Yli dj a n-j an d one obtains 
< s n+1 >= -n < a\ > - Yli' 1 < dj(Tn-j >= ~ n < °l >■ Thus dp = d\- J2T 3 < > A -J_1 dA = 
dX + J2T < s j+i > X~i~ l dX =< ipip* > since ipip* = X)o° s n X~ n . The normalization is built into 
this construction. It is clear that dp =< ipip* > dX cannot hold globally since < ipip* > should have 



poles at D + D* and the correct global statement follows from [25, [26], namely (dp/ < tpip* >) = dfl 
where dfl is the unique meromorphic differential with a double pole at 00 and zeros at D + D*. To 
see this note that dp/ < ipip* > will have zeros at the poles of ipip* (i.e. at D + D*) and a double 
pole at 00; this characterizes dfl. 



In summary we see that the quantity ip*ijj is seen to determine the Whitham hierarchy (2.4) and 



5 



the differentials dp, dE, dfl, etc. via 

dp dE dn 



dh (2.6) 



In the hyperelliptic case with say R(A) = IIi^^A — Aj) one can explicitly write down the 
differentials dp, Ml, etc. as dp = P(A)/V~R, P(A) = + fl^A*-* and dn = Q(A)/VR, Q(A) = 
6A9+1 + YfQbjA?-3. Now from (gl|) (#) <9 T dp = dxdfi = (dn/dp)d x dp. The moduli here are the 
branch points A,- which depend on the slow variables so writing out (J(t), multiplying by (A — Aj) 3 / 2 , 
and passing to limits A — > Aj yields the Whitham equations as equations for Riemann invariants Aj. 
in the form dxA^ = Vk(A\, • • • , A2 S +i)9xAfc where the characteristic velocities = (c^l/^P) I A=A fe 
have an elementary expression (cf. [0, ^]). One can also use the integrals of motion (~ /„ above) as 
moduli (see remarks later). 

3 ACTION AND PREPOTENTIAL 

The viewpoint of [^] for example is now to look at SW theory in terms of a renormalization 
group map {G,T,m, hi} — > {a,(/i), ^"(a,)} where G ~ gauge group, r ~ bare coupling constant, m ~ 
mass scale, hi ~ symmetry breaking vev's, aj(/i) ~ background "fields", and T ~ prepotential (with 
a,i = dT/ddi, etc.). This map is then decomposed as {G, r, m,h} — > {T,,dS m i n } — > {a,i(h), ^"(aj)} 
where S is a Riemann surface and dS m i n is a meromorphic one form on S such that ddS / dhi is holo- 
morphic. Since the genus of £ may be larger than rg* = rank(G) one must make some adjustments 



for which we refer to Roughly dS m i n ~ dXsw (modulo a few nonzero times - cf. [?, 22, 23]) 
and h ~ moduli as will be indicated. The first map to {£, dS m i n } involves quite general ideas and 
can be viewed via a spectral cover of some base space E(t) and an integrable system determined by 
a (matrix) Lax operator based on E(t) (generally E(t) ~ P 1 or an elliptic curve). Thus one looks 
at £ : det(t — L(z)) = which defines X as a ramified cover of E and the integrals of motion (or 
Casimirs) for the corresponding integrable system will be identified with moduli hi. This is in keeping 
with the spirit of averaging in Section 2. In this context dXsw ~ tdu where duo is (the) normalized 



holomorphic differential on E{t) and this will be clarified below (cf. [?, 13, [2^, |26[). Regarding the 
second map the action integrals cij = Sa % = £4. dXsw and af = Sb x = §g. dXsw for (Aj, Bi) a canon- 
ical homology basis will play a fundamental role. To assist in the determination of dXsw we recall 
that the monodromy for period integrals § c duj {dwi a basis of normalized holomorphic differentials) 
is often obtained from the physics and one will have PF equations for period integrals, etc. Thus it 
will be sufficient to have ddXsw I 'dhk = J2i a kjdujj so e.g. ddi/dhk = § A . J2 a kjdtOj = J2 a kj$ji = oiki 
and daf /dhk = § B , J2&kjdu!j = J2 a kjbji where (bij) is the period matrix for £ (cf. also 26] and 
remarks below). For hyperelliptic £ one often takes ddXsw I 'dvf- = dCj^ = (A k dA/y) ~ unnormalized 
holomorphic differentials (y = \/P, R = Ili 9+2 (A — Aj), k = 0, • • • ,g — 1) where Vk (1 < k < g) are 



suitable "flat" moduli (see e.g. pq, 24, E9[ for good illustrations and calculations). 



Now from [^] (cf. also ||) one works from a general action dS and builds in times T n from 
the beginning. For convenience the development here is in terms of hyperelliptic spectral curves £ 
which arise automatically from periodic Toda hierarchies (cf . |5|, || [l4|] ) ; physically the framework is 
appropriate for SU(n) N = 2 susy YM theories with matter for example (and others). Then the T n 
correspond to quasiclassical (slow) times and the prepotential J- will correspond to F = logarithm 
of the tau function for the dispersionless Toda hierarchy. For convenience the development here is in 
terms of hyperelliptic spectral curves £ which arise automatically from periodic Toda hierarchies (cf. 
[||, ^, Q). Thus let duji be canonical holomorphic differentials on £ 5 = T, g (h^) where i = 1, • • • ,g 
and k = ,K. The hk here are to represent the moduli (e.g. they could correspond to branch 



6 



points Afc for hyperelliptic situations). One has then canonical normalizations § A . dujj = <5jj and 
§ B , du)j = bij(hk)- From the point of view of Yang-Mills (YM) theory the spectral curve T, g is de- 
termined by the gauge group G with ~ (1/fe) < Tr(p k > for a suitable scalar field in the adjoint 
representation of G which breaks the original gauge symmetry down to U(l) r ° (cf. here 0, [b|, |3(| 
for details). From the point of view of integrability the (1/k) < Tr<p k > correspond to Hamilto- 
nians (integrals of motion of the integrable system) and one is back to Whitham dynamics. Now 
in [p3| one begins with a meromorphic differential dS on T, g defined (partly) by the requirement 
(•) (ddS/dhk) ^ J2i a kf — holomorphic differential, where everything depends on the hk- Exis- 
tence of solutions depends on the number K of moduli and one will assume here that K = g (one 
expects K > g generically) . We will anticipate hyperelliptic curves here (e.g. Toda situations) and 
assume S g to be given by A 2 = ni 9+2 (A — A») = R g {A) for example. Let us pick two points A± G H g 
and use A for a complex coordinate near A+ or A_ (A is correct here, not A). One defines dCt n as a 
solution to (•) satisfying in addition (1) dCl± n (n > 1) has a pole of order n + 1 at X± and no other 
poles with dQ± n (X) = ±[(A — \±)~ n ~ 1 +o(l)]d\ and (2) dtio has simple poles at both \± with residues 
±1. Such dtl n exist for K = g for example and one has then (••) (ddQ, n ) / dhk) — J2i cr fek c ^ t 'i- O n the 
other hand differentials d£l n , satisfying (1) and (2), but not necessarily (••), arise in Whitham theory 
with normalization § A _ d£l n = 0. Note that the conditions (1) and (2) plus normalization define the 
dQ n and then dfl n ~ dQ n + Xa cfduJi must hold with c" = § A , d£l n . One can write then re" = § B . dfl n 
and refer to [^] for the tau function of Toda theory. Then one could use as in Section 2 for 
averaging or work directly from an action expression (a nice exposition of Toda averaging is in Q) 
and end up with an action term dS(hk,Ti) = J2i a>idwi + Yl°?oo T n d£l n where the hk are moduli and 
the T n are slow variables. A priori a^, duJi, and d£l n depend on hk and T n and further analysis is 
mandatory. 

Thus to be more systematic in developing a quasiclassical tau function whose logarithm corre- 
sponds to the Seiberg-Witten (SW) prepotential one takes a solution dS(T n \hk) of (•) (presumed 
to exist without loss of generality since we have indicated how such creatures could arise from pe- 
riod integrals and monodromy) such that (*) {ddS / dT n ) ~ d£l n , which would make dS a kind 
of generating function for all solutions of (•). Thus write at the quasiclassical or averaged level 
a j = §A dS'i a f = §B dS and stipulate that (dcij/dT n ) = and (daf /dT n ) = § B , dfl n = K^{h). 

3 J 3 "' J 3 J 

Now expand dS as (f) dS ~ Y^oo u rn{T)d(l rn {h) with coefficients u m (T) independent of h. Putting 
(f) into (*) we get 

ddS x^fdum ~ ddVL m dh k \ , . 

= :" n " ^s^m)" 1 <3 ' 1) 

We will write now d n = d/dT n and then upon comparing d n dS ~ d£l n ~ ctti n — J2i cfduoi with 
flPD one gets (44) d n u m = 5 mn =>• u m (T) = T m along with the formula J2kd n hk (E-oo = 



J2kdnhk&ki = — cf. The last equation represents a version of the Whitham dynamics for hk (dki 
is defined as indicated). In particular the moduli hk are necessarily T dependent. We note that 
normally one thinks of Whitham equations as equations for d n £lk as in fl2.4| ) or as equations for 
branch points as indicated after fl2.6|) (here d n ~ d/dT n for slow variables T n ). The branch points 
are moduli in these hyperelliptic situations so in a general sense one can think of any nonlinear first 
order PDE in the moduli (e.g. hk) as Whitham equations. In the present situation one has in fact 
expressions (dd(l n /dhk) = Y^^ki^i with d n dS = dCl n — J2 c fdiOi whose combination leads to (44), 
i.e. to d n hk<Jki = — cf . At this point no restrictions on the number of moduli have been made (but 
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see below). We refer also to Q for another viewpoint on Whitham theory. 

Now (f) and (44) imply (using (44) dCl n ~ dQ, n + Yi cf(h)duji; cf(h) = § Ai dtl n ) that one 
has dS ~ Y^oo T m dfl m (h\T) ~ Y^oo (T m dQ m + T m Y{ c^duji) Integrating this along Aj cycles and 
using (44) again we get aj = § A . dS = Y^oo ^mC™ = Y^oo Tm § A . d(l m leading to the desired 
expression dS ~ Yi ctidtOi + Y^oo T m d^l m - Note that aj = aj(/i(T)) means that the a, seem to be T 
dependent but the T dependence via % actually cancels out since ddj/dT n = has been stipulated. 
This means that one can add the aj as additional independent variables to the set of slow times 
T n and write dS = dS(ai,T n ). Upon doing this we should add an equation for ddS/dcn to (*). 
This can not be postulated arbitrarily since dS is already defined and hence such equations must 
be derived. To do this we assume now explicitly K = g so that (aki) in (44) will be a square 
matrix. Then some calculation yields (dh^/ ddj\ T=c ) = {(Tjk)~ l {h) from which one can conclude that 
dS = Y°?oo T m dCl m = Yi cadiOi + Y^oo T n dfl n with ddS/dai = duji and d n dS = dCl n . We can also 
write now hk(ai,T n ) once the have been established as independent variables. Next one looks at 
the B-periods of dS in the form 

„ r / 9 00 \ 9 00 

af=f dS=f ]T Oidwi + ]T T n dn n = ]T Oibijih) + Tntf(h) (3.2) 

Then evaluating daf jdai = § B .{ddS /dcii) = § B . duii = bij(h) we see from the symmetry bu = bji 
that there is some function T such that a® = (dJ-/daj). This J-(ai,T n ) is called a logarithm of a 
quasiclassical tau function or a prepotential. The arguments are essentially the moduli of solutions 
to the original Toda like or low energy YM dynamical system with di(hk) parametrizing the moduli 
of the curve S g and T n moduli for coordinate systems in the vicinity of punctures. It turns out that 
the point T n = for all n is usually singular (note a^i = Y°^oo T m &ki vanishes when all T n = so 
does not exist). One can give explicit formulas for many quantities in terms of contour integrals 
which we omit here (cf. p3| ) but we do note that J- is a homogeneous function of degree two, namely 
^ = \ (Ei ai(dT/d ai ) + Y°?oo TndnT) . 



4 SYMPLECTIC GEOMETRY 



The study of TFT and LG theories in connection with integrable systems and Whitham dynamics 
has an extensive literature (we mention here only |?], [||, [IT], 12, 15, ^7], [H]]). There are many 



visible connections between F in the quasiclassical tau function r = exp(F) for KP/Toda theories 
and the prepotential T ' ■ Further, T and the LG (super)potential W for example are also related 
to the Hodge-Kahler special geometry in the N = 2 wonderland (cf. |l^, [H[]). Recently in || 
it was shown that the theory of the prepotential for N = 2 susy YM leads to WDVV (Witten- 
Dijkgraaf- Verlinde-Verlinde) equations as in TFT (cf. |l5|]). Also in |L7]] one uses duality theory 
as in SW theory to produce a Legendre transform relating x, J-, and the probability density |^| 2 
for the Schrodinger equation; this is suggested as a way to quantize geometry. Thus there are many 
directions and connections. 



In this section we discuss briefly some constructions from [26] which were in part spelled out 
more completely in ||. The paper |26| gives a particularly appropriate development of the theory 
of integrable systems in connection with N = 2 susy YM. The role of dXsw as an action-period 
creature is emphasized and for hyperelliptic Toda systems of the type arising in iV = 2 susy YM 
the theory seems quite complete. For systems of Calogero-Moser type with underlying elliptic curve 
there is much less information. In any event one outstanding feature of this work is to show how 
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much can be based on the BA functions ip and ip*. One begins with a universal configuration space 
which is a moduli space of Riemann surfaces with TV punctures P a (better thought of as marked 
points) and two Abelian integrals E and Q with poles of order at most n = (n a ) and m = (m a ) at 
the P a . One defines an n a -jet [z a ] nct of coordinates near a puncture P a to be an equivalence class 
of coordinates z a with z a = z' a if z' a = z a + 0(z™ a+1 ); this space of jets has dimension equal to n a . 
For [z] n a jet near Pi we define an Abelian integral Q as a pair (dQ, cq) where dQ is a meromorphic 
differential on the surface £ and Q = E- m CkZ k +c Q +R Q log(z) if dQ = d(^ m c k z k )+R^(dz/ z). By 
integrating dQ along paths one extends the Abelian integral Q holomorphically to a neighborhood 
of any point in £/{Pi, • • • , P/v}. The analytic continuation will depend in general on the path. 
Let us concentrate on iV = 1, with (n,m) = (n, 1) fixing the singularities for E and Q, and take 
dE ~ d(z~ n + 0{z))dz ~ dVL n with dQ ~ d(z~ l + 0{z)) ~ dp (necessarily Rf = = 0). For 
comparison purposes below we will take n = 3 later as well. The universal configuration space is 
defined as A4 g (n, 1) = {£, Pi, [z] n , E, Q} where £ is a smooth Riemann surface of genus g; there are 
then 3g — 3 parameters for the Riemann surface, n for the jet [z] or order n, and g + 1 each for P 
and Q, giving 5g + n parameters (note e.g. E has one parameter for ce and g for the holomorphic 
differentials which can be added without modifying the singular expansion. For local coordinates one 
chooses now 5g + n objects (P ~ Pi, z ~ Z\) 

T k = jRes P (z k QdE), (k = 1, • • • , n); r Al , E = <f dE; (4.1) 

T Bi,E = f dE; t a n = f dQ; t Bi ,q = f dQ; en = f QdE 

JBi JAi JBi JAi 

where i = 1, • • • , g and it is QdE = pdE which will play the role of d\sw here. 

Now let T> be the open set in A4 g (n, 1) where the zero divisors of dE and dQ do not intersect 
(i.e. the sets {7; dE^) = 0} and {7; dQ(^) = 0} do not intersect). It is proved that near each 
point in V the 5g + n functions ta u e-> T Bi,E, T Ai,Q, tb^q, a% have linearly independent differ- 
entials, and thus define a local holomorphic coordinate system. Further the joint level sets of these 
functions (omitting the a,) define a smooth g-dimensional foliation of T>, independent of the choices 
made in defining the coordinates themselves. Now Ai g = M g (n,l) can be taken as a base space 
for two fibrations M 9 and M; M 9 has fiber S 9 (T,) ~ J(S) = Jacobian variety (via the Abel map 
(71, • • • , 7 g ) — ► J2f=x Jp doJj) and M has fiber S (all over a point (£, P, [z] n ,E, Q) E -M g )- We con- 
sider only leaves A4 of the foliation of T> indicated above and look at the fibration M or J\f 9 over the 
base M. One wants to define a symplectic form l^m on N 9 . First, although E and Q are multivalued 
on the universal fibration their differentials are well defined on M. The idea here is that E and Q are 
well defined near Pi and their analytic continuations by different paths can only change by multiples 
of their residues or periods along closed cycles. But on a leaf of the foliation the ambiguities remain 
constant and disappear upon differentiation. Hence one has differentials 5E and 5Q on the fibrations 
which reduce to dE and dQ acting on vectors tangent to the fiber. One can trivialize the fibration 

with the variables a±, ■ ■ ■ , a g along the leaf A4 and e.g. E along the fiber. Then dQ coincides with 
(dQ / dE)dE where dE is viewed as one of the elements of the basis of one forms for M and the full 
differential is 5Q = dQ + J2i(dQ / 'da^dcii = dQ + S E Q (also one takes ce = 0). 

Now if one considers the full differential 5 (QdE) on J\f it is readily seen that it is well de- 
fined despite the multivaluedness of Q. In fact the partial derivatives d ai (QdE) along the base 
A4 are holomorphic since the singular parts of the differentials as well as the ambiguities are all 
fixed. In particular (d / dai)(QdE) = diOi where diOi is a basis of normalized holomorphic differentials 



§ A . dujj = 5ij with § B . dwj = bij. To see this note that it is implicit in (41) since by definition of 
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the Oj, (ddi/dcij) = 5ij = § A .{d{QdE)/daj) which implies (d(QdE) / daj) = du>j. Formally then one 
defines ujm on J\f 9 via 

(9 \ g 9 

Q{li)dE{ji) \ = £ SQfri) A d£( 7i ) = £ d ai A du>i (4.2) 

The first expression seems formally reasonable on N g and the last appears to be a calculation of 
the form um = 5(QdE) = Yf\dai A [d(QdE)/da,i] = J2fdai A du>i. Now go to KP for illustration 
and background. One can work with a (nonspecial) divisor (71, • • • , 7 S ) giving rise to quasiperiodic 
functions of t = (t n ) 1 < n < 00, of the form Ui^ n 1 < i < n — 2, 2 < n < 00, which arise 
as solutions of an integrable hierarchy. The BA function is defined as a meromorphic function 
away from P with simple poles at the 7,- (1 < i < g) and an essential singularity at P of the 
form ip(t,z) = exp tnZ~ n ) (1 + Ei° £,i(f) z )- There is a Lax operator (in general one for each 
puncture) L n = d n + J2o~ 2 u i.nd % with (d/dt n — L n )ip(t,z) = (d = d x , x = t\). Thus there is a 
map {£, Pz, 71, • • • ,7^} — > {ui tn (t)}. An explicit form for the BA function for KP is given in [26] 
but it involves normalizations 5ft § c dQ n = for any cycle C and is more complicated in appearance 
that our previous expression in Section 2 for example. Similarly the dual BA function tp* is defined 
as before (and denoted by ift in [p6|l). We note that for the Toda lattice one takes N = 2 punctures. 

An element in M 9 (n, 1) gives rise to a datum in a space Af 9 via H : (£, P, [z] n , E, Q, 71, • ■ • , 7 S ) — > 
(£, P,z,7i, • • ■ ,7 9 ) -» { ^(t)!™^ 2 } G Take a real leaf M (i.e. 5ft / c dE = 5ft jf c = for all 
cycles C on T) and write still dQ ~ <if2i = dp and <iE ~ dQ n (with real normalizations). We also 
take t\ ~ x. One wants now to express u>m in terms of forms on the space of functions {tti, n (t)}- 
First the Ui tTl (t) can be written in terms of the asymptotic BA coefficients £j and one knows that 
(d x ip/tp) = z' 1 + J2T h sz s and (d x ip*/ip*) = -z~ l + J2T K* s (<*• @ and arguments before Q 



where e.g. <9 = L + J2T a j^ ^ which implies e.g. (dtp/if)) = A + J2T a }^~^ f° r ^ ~ z ^ 1 )- I n 
any event the first n — 1 coefficients /i*, or crj are differential polynomials in the Ui tTl (initial 
data ^ s (t)\ x=0 = p> s (t2,-") determine £ s for s < n — 1). Writing H s =< h s > one gets now 
p = z~ x + H s z s (cf. @). We note that one uses < > x and < > xy averaging at various places in 



[26] but generically this should correspond to ergodic averaging. 



A result in now asserts that for M 9 the Jacobian bundle over a real leaf A4 of the moduli 
space A4 g (n, 1) the symplectic form uj_m can be written as 

^ = -i?es P Z—^—dp = nJ2 <Sh * A J 6 h n-s > (4-3) 



where the /i* are differential polynomials as above, and the differential forms 5h s and 5*h* are 
defined via Sh s = ^ 8 Ui , n Y,t{-d)Kdh s /duf n ) and 5% = E^O^SE^^ 



This is a fascinating result but the proof in [^] requires some embellishment. First one must 
come to terms with an expression (£6) 52£ = 5p(dE/dp) + (< ip*5L n ip > / < ip*ip >) where 
<5^n = Eo~ 2 fi u i,nd l (see below) and then it will follow from (|4.2|) (with 5p ~ 5Q) that (#•) ojm = 
— Ef(< "P*5L n ip > I < ip*ip >)(7s) A <5p(7 s ). The formula (44k) is asserted to come from |2j| but to 
see this one has to interpret 5v,i n as arising from edxUi^ n (X , T) in the first order term. Indeed we can 
look at ( [2.5D with the last two terms absent on the leaves of our foliation and written generically for 
L n ~ n n as (•*) < r{L l n d x L k - L\d x L n )^ >= d x Q k < i>*L\-p > -d x n n < ^L\^ >. Then for 
L n = L3 ~ and = L\ ~ p one has Li = 3, 5Li ~ d x L\ = 0, L\ = 1 and — < ip*d x L^ , p >= 
d x p < ip*L\ip > —d x Q,3 < ip*ip >. Written in terms of 5p ~ ed x p, SL^ ~ ed x L^ n etc. one obtains 
5fl 3 = 5p{< ip*L\il) > / < V> >) + (< 4>*SL 3 i; > I <^P*iP >) while from (§J) (with L\ ~ ^ x ) we 
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have (< i^*Llip > / < ip*tp >) = (dQ 3 /dp), or generically (dQ n /dp) = (< ijfLfy > / < >), 
and i&A follows. The formula (•♦) can also be used for a more general theorem in J26[ |. Thus 
(♦•) holds and writing 7 s (t) for the zeros of ifi away from P (corresponding to the fixed poles 
j s for t = 0) one can pick ti,---,t g (generically) as times for which the flows 7 s (i) are indepen- 
dent and use them as coordinates on S 9 (T,). These can be transferred to the system of coordinates 
/(7i)> • • • , /(7s)) for / an Abelian integral on S via (d/dti)f(j(t)) = Res^ t )[(d/dti)ip(t, z)/ip(t, z)]df. 
In the present case, writing (5 t tp/ip) = Xa^jVvVO^j an d using Res ls = —Resp one obtains 
u>M = Resp [(< tp*SL n ip > / < ift*ip >) A (<5; ^/V 7 ) lt=o] dp. It is not unnatural to see here an apparent 
change in the number of parameters. Now note that ip*5L n i/j = i[>* J2o~ 2 ^Ui,nd z tp = J2o~ 2 Sui^^d 7 " 1 ^ 
and an easy calculation gives (d^ip/ip) = z~i(l + J2T c jp zP )\ this implies (vp*(5L n ip)/ < ip*ip >)dp = 
— J2 C j SJ s z~ n ~ 1JrS dz for coefficients 5J S ~ linear combinations of the Sui iTl (recall (dp/ < tJj*tP >) = 
0(z~ 2 )). Averaging now gives (< ip*(5L n if) > / < ip*ip >)dp = < SJ S > z~ n ~ l+s dz so 

implies 5H S =< 5J S >= for 1 < s < n — 1. Hence ujm = ~J2i < $Jj+n > A dtj. Further from 
0, 0] one has (djt/)*/^*) = -X j -J2T = -z^+Ei° K,j z ' s and a calculation gives 5J j+n (t) = 

Yyy~ fiJsh>n-sji hence oj_m = —Resp(< 5t^* AdL n tp > / < ip*ip >)dp holds and an argument is sug- 
gested in p6[| to extend 5t to a full 5 so that the first formula in ((D|) is proved. For the second formula 
in Q one looks at blog^ = 6 (c(t i} i>2)+ f* o d^Zo*^*) = 5J2T (c 8 (*i, » > 2) + /* /*> s cfe) . 
This implies -Res P (< (dlogip*) A (if)*6L n ip) >)/ < ^ >) = o;^ = -Ei" 1 < ^Aj^ S*h* n _ s dx >. 
By definitions 5J S does not contain variations of derivatives of U{ so we can write 5J S = —n(dh s /du)5u 



and this gives the second identity in (4.3). We note here that there is a little interplay between 
5 ~ edx and X = ex in the last two lines. The examples in |2(| have some curious features but in 
any event one should keep in mind the conditions 5H S = when selecting functions. One notes that 
of course the symplectic structure involving 5h s , S*h* s , etc. does not include the time dynamics so 
KdV ~ L2 = d 2 + u for example and Bousinesq ~ L3 = <9 3 + ud + v. Adaptions and applications of 
the above fibration framework to SU(n), N = 2 susy YM theory are also given in 
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